Quantum discord amplification induced by quantum phase transition via a 
cavity-Bose-Einstein-condensate system 
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We propose a theoretical scheme to realize a sensitive amplification of quantum discord (QD) between two 
atomic qubits via a cavity-Bose-Einstein condensate (BEC) system which was used to firstly realize the Dicke 
quantum phase transition (QPT) [Nature 464, 1301 (2010)]. It is shown that the influence of the cavity-BEC 
system upon the two qubits is equivalent to a phase decoherence environment. It is found that QPT in the 
cavity-BEC system is the physical mechanism of the sensitive QD amplification. 
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Quantum discord (QD) 01 01 is considered to be a more 
general resource than quantum entanglement in quantum in- 
formation processing Jj-^]. A nonzero QD in some separa- 
ble states is responsible for the quantum computational effi- 
ciency of deterministic quantum computation with one pure 
qubitHHi and also has been considered as an useful re- 
source in quantum locking ||5|] and quantum state discrimi- 
nation IslZI]. On the other hand, any realistic quantum sys- 
tems interact inevitably with their surrounding environments, 
which introduce quantum noise into the systems. As an use- 
ful resource, it is interesting that QD can be amplified by the 
quantum noise. We find the QD does be amplified for two 
non-interacting qubits immersed in a common phase decoher- 
ence environment II 1 Ofl - Especially, when the two qubits are 
identical, the phase decoherence can induce a stable amplifi- 
cation of the initially-prepared QD for certain X-type states. 
In this paper, we propose a scheme to realize the controllable 
QD amplification of two atomic qubits by making use of an ar- 
tificial phase decoherence environment consisting of a cavity- 
Bose-Einstein condensate (BEC) system. 

The Dicke model lUlllT^ll describes a larger number of two- 
level atoms interacting with a single cavity field mode. As in- 
creasing atom-filed coupling, the model predicts a QPT lfl3ll 
from the normal phase, which the atoms are in the ground state 
associated with vacuum field state, to the super-radiant phase, 
which both the atoms and field have collective excitations. 
Recently, a cavity-Bose-Einstein condensate (BEC) system is 
employed to firstly realize the Dicke quantum phase transi- 
tion (QPT) experimentally and to explore symmetry breaking 
at the Dicke QPT 13. Meanwhile, the QPT system usually 
displays ultra-sensitivity in its dynamical evolution near the 
quantum critical point JH5l - ll9ll . which has been confirmed by 
a NMR experiment Il20ll . The purpose of this paper is to show 
that the QD of two initially correlated atomic qubits can be 
sensitively amplified via the cavity-BEC system near the crit- 
ical point. We show that the cavity-BEC system can form an 
artificial phase decoherence environment for the two atomic 



qubits, and the QD of the two atomic qubits can be amplified 
by adjusting the QPT parameter of the cavity-BEC system. 

The physical system under our consideration is shown in 
Fig. HI A BEC with N identical two-level &7 Rb atoms is con- 
fined in a ultrahigh-finesse optical cavity. The atoms inter- 
act with a single cavity model of frequency u> c and a trans- 
verse pump field of frequency to p . We consider a situation 
that the frequency <d c and a> p are detuned far from the atomic 
resonance frequency cja so that the dutunings far exceed the 
rate of atomic spontaneous emission, the atoms only scatter 
photons either along or transverse to the cavity axis. Before 
the pump field turns on, atoms in the BEC are supposed to 
be in the zero-momentum state \p x ,p z ) = |0, 0). As soon 
as one turns on the pump field, via the photon scattering of 
the pump and cavity fields, some atoms are excited into the 
momentum sates \p x ,p z ) = | ± k,±k) = ~Zi VlU2=± \ \v\k,V2k) 
(hereafter, we take % = 1.) due to the conservation of mo- 
mentum, where k is the wave-vector, which is approximated 
to be equal on the cavity and pump fields. The two momen- 
tum states with |0, 0) and | ± k, +k) are regarded as two levels 
of the &7 Rb atom with energy separation o)q = k 2 /m, where 
m is the mass of &1 Rb atom. Defined the collective operators 

I = £; | ± k, ±k){+k, +k\, J + = Jl = E ; I ± k, ±k)(0, 0| with 
the index i labelling the atoms, the cavity- BEC is described 
by the Dick model lfl4l|2lll 



Hi — (jja^a + ojqJ z -\ — (a + 5^) (7+ + /_), 



(i) 
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where a 1 (a) is the creation (annihilation) operator of the cavity 
field, oj = -A c + U[)N/2 is effective frequency of the cav- 
ity field including the frequency shift induced by the BEC 
under the frequency of pump field rotating frame, where 
A c = oj p - <d c is the dutuning between pump field and cav- 
ity field and Uo = #q/A is the frequency shift of a single atom 
with maximally cavity field coupling strength go and detuning 
A = oj p - oja- A = '\fNgoQ.p/2A is the coupling strength in- 
duced by the cavity field and pump field with Q p denoting the 
maximum pump Rabi frequency which can be adjusted by the 
pump power. 

We consider such a situation that the two atomic qubits pass 
through the cavity at the same time and interact with the single 
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FIG. 1: (Color online) Schematic of our physical system: Two 
atomic qubits A and B with energy separation coa and a> B are, re- 
spectively, injected into the cavity in which a atomic BEC couples to 
a single cavity field and a transverse pump field. 



cavity field, the expression of the Hamiltonian reads as 



H 2 



^a+^^ + ^^gAtfa* + gB aW_ + H.c) , (2) 



where & z <B> = \e) (e\ - \g) (g\ is pauli operator with \e) and \g) 
being the excited and ground states. o% <B) (cr A( ^) is the rais- 
ing operator(lowering operator). gA(B) is the coupling strength 
between the atomic qubit A(B) and the cavity field, u>a(b) is 
the energy separation. Here we have made a rotating wave 
approximation. If the atomic qubit is far-off-resonant with the 
cavity field satisfying the detuning A A(B) = o>mb) - <w is much 
large than the corresponding coupling coupling stre ngth g A (B), 
one can use the Frohlich-Nakajima transformation II221 12311 to 
make the Hamiltonian in Eq. (f2j) become the following ex- 
pression 



H 2 = \to A &^ + \tu b <t b z + (cj + 6 A cr A + £ B <xf)a t a, 



(3) 



where aj' A(B) = uj A(B) + 6 A (B) with 5 A (B) = g 2 A{B) I^A(B) being 
the frequency shift induced by the scattering between cavity 
field and atomic qubit A(B). Then the effective Hamiltonian 
describing the two atomic qubits passing through the cavity- 
BEC system is 



H e ff = \jW A o^ + x« B a-f + {Sao^ + 6 B <r B ) a'^a + H\ . 



(4) 



Now we consider the dynamics of the two atomic qubits 
passing through the cavity-BEC system. We assume the two 
atomic qubits are initially prepared in a class of state with 
maximally mixed marginals (p A(B) = 7 A(B) /2) described by 
the three-parameter X-type density matrix p s (Q) = 1 /4(I AB + 
2; =] ci&f <S> of ), where I AB is the identity operator in the 
Hilbert space of the two atomic qubits, i = 1,2,3 mean x, y, z 
correspondingly, and c, (0 < |c,j < 1) are real numbers satisfy- 
ing the unit trace and positivity conditions of the density op- 
erator p,(0). The cavity-BEC system is initially in the ground 
state |G) of the Hamiltonian in Eq. ([T). The dynamic evolu- 
tion of the total system is controlled by the Hamiltonian in Eq. 



©. The density operator of the total system at time t is writ- 
ten as p T (t) = U(p s (Q) ® \G)(G\)W with U = e^ff'. After 
tracing the degree of freedom of the cavity-BEC system, we 
obtain the reduced density of the two atomic qubits 



l+c 3 p(t)Di(t)^ 

1 - c 3 v(t)D 2 (t) 

v*(t)D* 2 (t) l-c 3 

fi*(t)D*.(t) 1 + c 3 



(5) 



where we have introduced the following parameters 

H(t) = (ci - c 2 )e- i(£ ^ +a, * )( , v(r) = (ci + C2)e- i( -^- b) '" )t , 



with 



iH.,t-iH„t 



e -'"si'\G), D 2 (t) = {G\e iH <s'e- iH °''\G), (6) 



-6\a'a + H\,6\ — &a + <5b 
-8 2 a i a + H\,S 2 = 5 A - 5 B - (7) 



We consider the situation that two atomic qubits pass 
through the cavity field region in a very short time of satis- 
fying the conditions S\t <K 1 and 6 2 t <k 1. In fact, according 
to Ref. [ 14], the waist of the cavity field is 25 /urn, the effective 
frequency shift 6a, 6b are about 100 Hz, above conditions are 
well satisfied if injected velocity of the atomic qubits meets 
v » 10~ 3 m/s. By the short time approximation, the factors 
l-Di(0l> IAz(0l can be derived as 



\Dx if) I 



exp (-2ySy) , \D 2 (t) | = exp (~2y6 2 2 t 2 ) . 



(8) 



where the decay factor y — (ja^aj J - (a^aj is the cavity 

photon number fluctuation (PNF) in the ground state |G) lfl9ll . 

From Eq. (5) we can see that the cavity-BEC system only 
affects off-diagonal elements of the density for the two atomic 
qubits, hence it is equivalently a phase decoherence environ- 
ment for the two atomic qubits. That is, the cavity-BEC sys- 
tem constitutes an artificial phase decoherence environment 
of the two qubits. The QPT parameter of the cavity-BEC sys- 
tem A is a controllable parameter of the artificial environment. 
It's worth noting that when the effective frequency shift 6a, 6b 
are equal, i.e, 6 2 — 0, a decoherence free space in the basis 
i\ge), \eg)} appears. 

In order to obtain the detailed form of the PNF y, in the 
following we give the ground state \G) according to the Ref. 
lfl2Tl . Utilizing the Holstein-Primakoff transformation ll24ll 
J + = c T y/2 j - c'c, J- = y2 j - c'cc, J z = c'c - j, where j = 
N/2 , the Hamiltonian of the Eq. (Q]) is further reduce to 



Hi — wa'fl + o»oc'c + A.(a + a') 



2j 



+ H.c. 



(9) 



When the coupling strength A is smaller than the critical cou- 
pling strength A c = ywwo/2, i.e., A < A c , the system is in the 
normal phase where the BEC and the cavity field have low ex- 
citations. While when the coupling strength is larger than the 
critical strength, i.e., A > A c , the system is in the super-radiant 
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phase where both the BEC and the cavity field have collective 
excitations in the order of the atom number N. 

In the normal phase at the thermodynamic limit j — > oo, 
neglecting terms with j in the denominator, Hamiltonian (0 
becomes 



Hi = coa ' a + ojqc + A [a + $\ [c + c 1 j , 



(10) 



where we omit the constant term. The Hamiltonian in Eq. ( TTOb 
can be diagonalized as 



H\ — E-d[d\ + s+dLdq. 
by the Bogoliubov transformation 

a 1 = f\d\ + hd\ + fd[ + f 4 d 2 , 
c f = h\d\ + hjdi + h^d 2 + h^di, 



(11) 



(12) 



where the eigenfrequencies £_ and e + of the cavity-BEC sys- 
tem have the following expression 



* = I 

e± 2 



+ Wq ± -J (u)q - ai 2 ^ + 16A 2 ojcjo 



(13) 



The coefficients of Bogoliubov transformation about the cav- 
ity field in the normal phase are 

1 cos0 1 sin0 

/U = --=L(W±C_), /3,4= =-=L(w±C + ), (14) 

2 ye_o> 2 ye+w 

where the mixing angle is given by tan 2<p = 4 ^J^T ■ 

In the supper-radiant phase, we displace the bosonic modes 
a ' — > a" + yfa, c + — > - -y^S with \/a and -y^S describ- 
ing the macroscopic mean fields above A c in the order of 
O(j). Neglecting terms with j in the denominator and tak- 

ing = 77 4{iX-f), ^ = V7(W) with^ = | , the 
Hamiltonian Eq. © is reduced to the following form 

,\2 



Hi — coa ' a + cjqc ' c + r] (c ' + c ) 
+~A(a + a 1 ') (c 1 + cm , 

where the parameters (bo, A and 77 are given by 



(15) 



— (1+40. A = /^4 : -. , 

2£ v 97 \l+4 

w (1- 4-) + 40 



(16) 



= s_d l d l + s + d 2 d 2 



84" (i+ 4-) 

The Hamiltonian in Eq. (TT~5T > also can be diagonalized as 

(17) 

by the Bogoliubov transformation 

a Y = f[d\' +f' 2 di + f' A S 2 , (18) 

c' = h x d\ +h 2 di+hi t d\ + h 4 d 2 , 




FIG. 2: (Color online) The logarithm to base 10 of the PNF 7 changes 
with the coupling strength A. Related parameters N = 10 5 , a> = 
0.05MHz, <±> = 20MHz correspond the experimental parameters in 
Ref. (H. 



where the eigenfrequencies s_ and s + read as 



2 ^ w o 



e2 



•4w 2 t 



(19) 



The coefficients of Bogoliubov transformation about the cav- 
ity field in the super-radiant phase are 



A: 



1 cos 
2 



(w ± s'_) , f X4 



1 sin</> 



(w±e' + ). (20) 



2o>o>o£ 2 



where is the mixing angle defined by tan 2<p 

The PNF can be given in the normal phase with ground 
state 10,0)^42 and the super-radiant phase with ground state 
[0, 0) d > j , respectively, as the following forms 



7 = 



2ft si + 2fift + (/1/4 + fihr . a < A,, 



(21) 



U+/ 2 ) 2 + (/3+/J 



A > Ac. 



Compared with the case of the normal phase, the displace- 
ment a due to collective excitation appears in the super- 
radiant phase. Figure|2]shows the PNF y will experience dras- 
tic change near the critical coupling point A c /o»o = 10. The 
closer the coupling strength A near the critical coupling point, 
the larger the PNF y. This inspires us to control the coherence 
decay rate of the two atomic qubits by adjusting the pumping 
power to change the coupling strength in the region near the 
critical coupling. 

In the following we consider the QD amplification of the 
two atomic qubits induced by the QPT of the cavity-BEC 
system. The QD fl]] is defined as the difference between 
the total correlation and the classical correlation with the ex- 
pression T>(p AB ) = l(p A : p B ) - c{p AB ) with p A , p B , and 
p B being the reduced density operators for subsystems A and 
B, and the total density operator, respectively. The total 
correlation in the state p AB is measured by quantum mutual 
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FIG. 3: (Color online) The QD amplification rate as a function of the 
coupling strength A and the initial parameter c\ . Other parameters are 
set as c 3 = Ci/2, c 2 = 0, co = 0.05 MHz, tu = 20 MHz, t f = l/coo, 
Si = 0.001 wo, d 2 = 0, and N = 10 5 . 

information I [p A : p B ) = S (p A ) + S (p B ) - S [p AB ] with 
S (p) — -Tr(p log p) being the von Neumann entropy. The 
classical correlation between the two subsystems A and B 
is given by C(p AB ) = S(p A ) - min { ^, p k S (p£)] where 
Pk = Tr AB [(/ A ®^)p AB (/ A ®^)] denotes the probability relat- 
ing to the outcome k, and I A denotes the identity operator for 
the subsystem A with \P B ) being a set of projects performed 
locally on the subsystem B. 

The mutual information of the state given in Eq. © is 
derived as l(p A : p B ) = 2 + 2^ =1 /l/ log /I/, where Ai t 2 = 

1(1 + c 3 ± ^3,4 = |(1 - c 3 + |y(0£»2(0l) are four 

eigenvalues of p s (t). And the classical correlation can be ob- 
tained as JTSIH C(p s (t)) = | i^Y^ log 2 [1 + {-l) n x] with 

*(r) = max[|c 3 |,(IM0£>i(f)r+|v(f)£>2(f)l)/2]. Therefore, the 
QD can be written as 

4 

Dip s {t)) = 2 + ^,log 2 ^-C(p. s (f))- (22) 

i=l 

The QD can be amplified for some initial states such as the 
state parameters being set as C2 = 0, < c\ — 2c3 < 2/3 



when the qubits are in the phase decoherence environment 
ifioll . For the present cavity-BEC environment, let the two 
atomic qubits enter the cavity at time t — and leave the cav- 
ity at time tf. Then we can define the QD amplification rate 
as T = D(tf)/D{0). In Figure|3]we have plotted the amplifi- 
cation rate T with respect to the coupling strength A and the 
initial state parameter c\ when C3 = ci/2, C2 = 0, wo = 0.05 
MHz, cd = 20 MHz, AJuq = 10, t f = l/co , 6 X = 0.001 cj , 
6 2 = 0, and N = 10 5 . Fi gure [3] indicates that the initial QD 
can be amplified by the use of the cavity-BEC system through 
changing the QPT parameter A. Specially, the QD amplifica- 
tion rate sensitively increases at the QPT point of the cavity- 
BEC system A = A c . In this sense, the sensitive QD amplifi- 
cation can be understood as a quantum phenomenon induced 
by the QPT of the cavity-BEC system. It should be pointed 
out that one can control the QPT parameter A by changing 
the Rabi frequency of the pump field Q p due to the relation 
A = y/Ng n p /2A. 

In conclusion, we have proposed a scheme to realize the 
sensitive QD amplification of two atomic qubits via the cavity- 
BEC system through by changing the QPT parameter of the 
the cavity-BEC system, and and revealed the QPT mecha- 
nism of the sensitive QD amplification. We have indicated 
that the cavity-BEC system is equivalent to a phase decoher- 
ence environment for the two atomic qubits. Hence, it pro- 
vides an artificial and controllable phase decoherence envi- 
ronment for quantum information processing. It should men- 
tioned that the present scheme should be within the reach of 
present-day techniques since the cavity-BEC system used in 
the scheme has been well established in recent experiments of 
observing the Dicke QPT [14]. The experimental realization 
of the scheme proposed in the present paper deserves further 
investigation. 
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